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ABSTRACT _ 

Ammonia is a common pollutant in water as the result of runoff in 
agricultural areas where it is applied as fertilizer. It must be monitored 
regulary for safety purpose. The current testing technique does not allow on¬ 
site measurement as the equipment are bulky, the measurement process is 
time-consuming and tedious with additional mixing reagents. In this study, 
the presence of ammonia in distilled water is estimated from open-ended 
coaxial probe in the range of 200 MHz to 14000 MHz. Experimental results 
were obtained from two set of samples with seven different ammonia 
concentration each. The measurements are repeated thrice hence producing 
forty-two data sets with 550 points. Both curve fitting and multiple 
regression analysis were considered to perform valid ammonia concentration 
projection. Validation based on 5-fold and 10-fold cross validation suggested 
the feasibility of the technique as it presents root mean square error (RMSE) 
which is less than 0.02 in the ammonia prediction. Detection method based 
on open-ended probe would be convenient, simple and accurate method for 
in-situ determination of ammonia concentration. 
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1. INTRODUCTION 

Nitrogen plays a vital role in the formation of amino acids which is required by the plant for 
survival. However when it exists in molecular form it cannot be absorb by the plant directly, so it must be 
converted. As the results of conversion, ammonia (NH3) is produced [1]. Ammonia is a colorless compound 
consists of hydrogen and nitrogen that is highly soluble in water. It is a common pollutant in sewage, landfill 
leachate, palm mill effluent and industrial effluent [2, 3]. Ammonia caused eutrophication, acidification of 
soil, disturb the water treatment plant operation, and animal mortality [4]. While, for human ammonia may 
cause symptoms such as inhalation and indigestion problems [3, 5, 6]. 

Monitoring of water quality is an important aspect in the process of determining the pollution level 
of water. Through water quality monitoring, the researchers can predict the condition of environment caused 
by anthropogenic or natural activities. The data analysis is able to determine the suitability of remediation 
effort to ensure environmental standards are met. 

In Malaysia, the monitoring program on the safety of drinking water is handled by Department of 
Environment, Ministry of Natural Resources and Environment. The testing procedures for the water is well 
documented in the national standard [7]. The monitoring process starting by sampling, storage, handling and 
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analysis in the laboratory. The samples taken must be kept in sealed glass or plastic bottle. Once they are 
collected, they must be placed in sturdy boxes, packed with ice and kept in a dark environment at constant 
temperature in between 4-10°C. The sample box shall reach the laboratory within 24 hours seal with 
appropriate pre-treatment and great care [7]. All these steps are cumbersome and the integrity of the samples 
might not be preserved. 

Laboratory testing method such as indophenol reaction is only applicable when the ammonia 
concentrations is in between 0.025 to 3 mg/L. The testing results relies on the observation (changes of colour 
of the mixture through the titration), the uncertainty of the testing results is high as the colour representation 
is valid for a range of ammonia concentration [8]. Reagent test through spectrophotometer requires additional 
substance to be added into water and the ammonia concentration is judged based on penetration of light 
through the sample [9]. This equipment is expensive and must be calibrated regulary to ensure 
the reliability [10]. Ammonia-selective electrode involve tedious procedure and standard operating procedure 
and less sensitive [11]. Overall, the current testing procedures of ammonia is time-consuming, expensive, 
tedious and cumbersome. Due to this reason, alternative approach which enable the measurements to be 
performed at the location of interest is desirable. This paper intend to explore the feasibility of open-ended 
coaxial probe for ammonia concentration determination. 

Open-ended coaxial probe is one of the dielectric measurement techniques [12]. It is used to 
determine the relative permittivity of dielectric materials. Relative permittivity is a complex parameter which 
can be represented by (1). The real part of (1) is known as dielectric constant, it measures the energy level 
stored in a material when it is exposed to external electric field. While imaginary part of (1) is known as loss 
factor. This parameter represent how dissipative or lossy a material is in the presence of an external electric 
field [13]. Since every material is unique, so they will have individual set of electrical characteristic. 
Based on this, it is believed that water with different ammonia concentration will possess different relative 
permittivity. 

£ = €r-j€i ( 1 ) 

In this context, samples with different ammonia concentration are prepared. Their loss factor (LE) 
are measured based on open-ended coaxial probe. Based on the measurement data, curve fitting and multiple 
regression analysis are conducted to produce mathematical model which can predict the ammonia 
concentration accurately. The following sections will describe in detail the samples preparation, dielectric 
measurement setup and the formation of mathematical model. The efficiency of the model is confimed 
through the cross validation. 


2. SAMPLE PREPARATION 

Ammonium chloride has been chosen for sample preparation as it is a cheaper solution and it can be 
obtained easily compared to the pure ammonia liquid. Ammonium chloride is combination of ammonia with 
hydrochloric acid as shown by the chemical formula (2). Samples with ammonia concentration of 0 g/L, 0.1 
g/L, 0.2 g/L, 0.4 g/L, 0.6 g/L, and 0.8 g/L are prepared. These sample are acquired by diluting a stock 
solution which contains 1 g/L of ammonia concentration. In order to prepare this stock solution, it is 
important to identify the molar mass of the elements. The molar mass of the ammonia, hydrochloric acid, and 
ammonium chloride are 17.031 g/mol, 36.46 g/mol and 53.491 g/mol respectively [14]. 

NHs HCL=NH4CL (2) 

The number of mole of ammonia can be obtained by dividing the mass of ammonia required with 
the molar mass of ammonia: 


Number of mole of ammonia 


Mass of ammonia required 
Molar mass of Ammonia 


^9 

limig/l 


0.0587 mole 


The mass of ammonium chloride required in order to obtain stock sample with ammonia 
concentration of 1 g/L can be obtained by multiplying the number of mole of ammonia with the molar mass 
of ammonium chloride. So, 3.1408 g of ammonium chloride will be mixed into 1 litter of distilled water for 
stock solution preparation. 


Mass of ammonia = number of mole of ammonia x molar mass of ammonium chloride 
= 0.0587 mole x 53.491^ = 3.14085 
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After the stock solution has been prepared, it is further diluted in order to prepare more samples with 
lower concentration. Seven samples with concentration of 0 g/L, 0.1 g/L, 0.2 g/L, 0.4 g/L, 0.6 g/L, 0.8 g/L 
and 1.0 g/L have been prepared based on (3), where Mi represent the concentration of the stock solution, 
Vi represent the volume of stock solution that should be use to be diluted by V 2 , M 2 is the desire 
concentration and V 2 is the volume of the desire solution. Table 1 list down the required volume that should 
be mixed into 150 ml of distilled water in order to produce the desire samples with specific concentration. 
Two set of samples are prepared to ensure the measurement uncertainty from measuring scale and cylinder 
will not affect the concentration of the samples. 

= M2V2 (3) 


Table 1. Volume from stock solution required for samples prepa ration 
Concentration of ammonia, g/L Volume from stock solution, ml 


0.1 

15 

0.2 

30 

0.4 

60 

0.6 

90 

0.8 

120 


3. OPEN-ENDED COAXIAL PROBE FOR DIELECTRIC MEASUREMENT 

The open-ended coaxial probe has been chosen for LF measurement as it is more suitable and 
convenient for liquid samples compared to the other techniques [15]. This measurement setup consists of an 
open-ended coaxial probe, it is connected to vector network analyser (VNA). The VNA is connected to a 
computer embedded with specific software for the display of LF results, as shown in Figure 1. 

Before the measurement begin, the probe will go through open, short and load calibration process. 
This process is to ensure the uncertainty due to the cable, and connectors will be eliminated and the reading 
from the VNA is purely relies on the sample under test (SUT). During the measurement, the probe will be 
immersed into the SUT as shown in Figure 2. 



Figure 1. The loss factor measurement setup based Figure 2. The probe is emerged into the SUT 
on open-ended coaxial probe 


It is important to ensure there is no air bubbles exist on the probe as this will reduce the 
measurement accuracy. Two set of samples will be measured thrice separately in order to ensure the 
reliability of the loss factor results. The loss factor results will be discussed further in the result and 
discussion. 


4. PREDICTION MODEL BASED ON CURVE FITTING AND MULTIPLE REGRESSION 

The measurements were perfomed with six repetitions where each repetition consists of seven data 
sets at different concentration (0 g/L, 0.1 g/L, 0.2 g/L, 0.4 g/L, 0.6 g/L, 0.8 g/L and 1.0 g/L). Thus, in total 
there are forty-two data sets representing multiple LF of samples with different ammonia concentration. Each 
of the LF data set consists of 550 discrete data ranging from 200 MHz to 14000 MHz as shown in Figure 3. 
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Frequency, GHz 


Eigure 3. Loss factor of samples with different ammonia concentration 


Eor the purpose of deriving a mathematical model that best describe the characteristics of these data, 
we have attempted curve fitting. We have assumed that by applying the outcome of the curve fit, LE can be 
feasibly predicted at any frequency points. It is found that the simple graph models such as linear, 
polynomial, exponential, logarithm and power may not fit the data well. Hence, we have proposed a non¬ 
linear logarithmic curve function in (4) based on our best knowledge of graph theory and pattern heuristics. 
Erom the equation, f represents the frequency in MHz and LE represent the loss factor. At the end of the 
fitting process, forty-two sets of “b” (bl, b2, and b3) or coefficient values shall be obtained. Based on these 
coefficient values, we have applied a simple multiple linear regression analysis to predict ammonia sample 
concentration (C). The generalized equation for multiple linear regression are shown in (5). 


LF = ln(/) + b 2 f + b^ 

(4) 

C = biQi + b2a2 + + d 

(5) 


5. RESULTS AND ANALYSIS 

5.1. Loss factor of water with different ammonia concentration 

Eigure 3 shows an example of the LE results among the six repetitions. In general, it is observed that 
samples with different ammonia concentration have different LE. The trend of the loss factor varies 
proportionally with the ammonia concentration. Nevertheless, the pattern is significantly different especially 
at frequency lower than 4 GHz. As the frequency increases, the variation of the LE become lesser even 
though there are changes in the ammonia concentration. It is interesting to found that sample with 0 g/L and 
0.1 g/L ammonia concentration has similar loss factor in between 8.5 GHz to 10.5 GHz. This condition 
indicates that these two samples are similar in terms of electrical characteristics. Specifically, we also found 
that LE is insensitive on the ammonia concentration variation at higher frequency. 

5.2. Prediction results based on curve fitting and multiple regression 

Eorty-two data sets had been modelled by non-linear logarithmic curve function (4). The values of 
fitting variables to is partially summarised in Table 2. It is observed from Table 2 that even for sample 
with same concentration, different values of were obtained. Such a difference is mainly due to the 
uncertainty in dielectric measurement and sample preparation, hence different measured LE values are 
obtained. We have showcase an example of curve fitting in Eigure 4. Erom the figure, we have found that in 
general our proposed curve-fitting function as in (4) is able to describe our experimental dataset well, albeit 
with some inaccuracy at the lower frequency range. 

Next, multiple regression analysis is performed by taking all the coefficient values of bi, b 2 , and bs 
in Table 2 as predictor for ammonia concentration. According to the generalized equation for multiple linear 
regression (5), the outcome of four coefficients values of aj, a 2 , as and intercept d are summarized in Table 3. 
Based on the outcome of correlation analysis, we have discovered that the multiple linear regression model as 
in (5) is a good predictor of ammonia concentration with a correlation coefficient (R^) of 0.998. A closer 
examination on the relationship of each parameters is as summarized in Eigure 5. Erom the figure, we have 
found that the pairing relationship among the variables are strong correlated as they are generally above 0.8 
for values. 
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Table 2. Partial results of fitting variables with corresponding concentration 


Data Sets 

bi 

62 

63 

Concentration (g/L) 

1 

-0.1939 

0.0031 

4.67 

0.1 

2 

-1.5434 

0.0033 

16.1922 

0.2 

3 

-3.3398 

0.0028 

33.0456 

0.4 

4 

-6.1733 

0.0031 

56.648 

0.6 

5 

-10.411 

0.0036 

89.7221 

0.8 

6 

-13.3015 

0.004 

112.355 

1 

7 

4.2951 

0.0018 

-26.9456 

0 

8 

2.9546 

0.0015 

-15.6748 

0.1 

9 

0.6524 

0.002 

2.0435 

0.2 

10 

-3.0014 

0.0025 

30.9112 

0.4 



Figure 4. Comparison between actual and fitted data at ammonia concentration 0.6 g/L 


Table 3. Multiple regression results 
Variables _ oi _ 02 _ 03 _ d 

Values 0.342 32.525 0.050 -0.118 

Observation: 42 
R^: 0.998 

Residual Standard Error: 0.018 


To further examine the merits of our proposal, we have performed ^-fold cross validation. In this 
research, we have examined the validity of our prediction model using 5-fold, 10-fold and 10-fold with two 
repetitions. The outcome of our cross validation are summarised in Table 4. In an overview, we have found 
that Root Mean Square Error (RMSE) is generally below 0.02, which indicates good fitting of dataset and the 
validity of our prediction model. The results also implicitly indicate that while multiple coefficient values 
are obtained from curve fittings, such a result diversity does not impact the final prediction accuracy of our 
model. In short, our proposed prediction model is robust in predicting the concentration level of ammonia 
samples based on empirical dataset. 

We wish to highlight that our model is only valid for ammonia concentration prediction from 0 to 1 
g/L with the frequency range from 200 MHz to 14000 MHz when the ammonium chloride is used. While we 
have not attempted with concentration above 1 g/L, we do believe that extrapolation of this model is possible. 
Secondly, while our approach has the potential benefit of shortening the field experimental process, 
our model is only limited to one type of contaminant and may not be valid if the sample contains other 
contaminants. Given our controlled experimental environment, our prediction model need to be further 
examined in the case of actual field experimental samples that well possess other types of contaminants. 
All these are potential areas for future research that we shall investigate further. 
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Eigure 5. Parameters pairing plot 


Table 4. Multiple regression results 


Validation type 

RMSE 

R" 

5-fold 

0.0178 

0.9980 

10 -fold 

0.0172 

0.9987 

10 -fold, 2 repetition 

0.0180 

0.9988 


6. CONCLUSION 

In this work, a predictive model for ammonia concentration estimation in water has been presented. 
This model is presented using curve fitting and multiple regression analysis of forty-two experimental LE 
data sets for ammonia concentration ranging from 0 g/L to 1 g/L. Erom our results, our proposed non-linear 
logarithmic equation can be modelled successfully. Eor the multiple regression model, the outcome of our 
5-fold and 10-fold cross validation results (RMSE) of less than 0.02 has successfully indicate the validity of 
our model. We reckoned that our prediction model can potentially reduce the ammonia testing duration 
compared to the conventional testing method conducted in laboratory, where dielectric measurement can be 
performed on-site using portable open-ended coaxial probe. Eor future work, we shall explore additional LE 
data for samples at higher ammonia concentration as well as predictive modelling in the presence of 
multiple contaminants. 
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